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Abstract 

Kaons and protons carry large parts of two conserved quantities, strangeness and baryon number. 
It is argued that their correlation and thus also fluctuations are sensitive to conditions prevailing 
at the anticipated parton-hadron phase boundary. Fluctuations of the (K + +K~)/(p + p) and K + /p 
ratios have been measured for the first time by NA49 in central Pb+Pb collisions at 5 SPS energies 
between y/s~^= 6.3 GeV and 17.3 GeV. Both ratios exhibit a change of sign in <7d y n, a measure of 
non-statistical fluctuations, around y/s NN = 8 GeV. Below this energy, Odyn is positive, indicating 
higher fluctuation compared to a mixed event background sample, while for higher energies, Odyn 
is negative, indicating correlated emission of kaons and protons. The results are compared to 
hadronic transport model calculations which fail to reproduce the energy dependence. 
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Heavy-ion collisions serve as the laboratory to study hadronic matter under extreme 
energy density and temperature conditions. They surpass critical values of energy den- 
sity where lattice QCD calculations predict a phase transition from hadronic to deconfmed 
matter (the "quark-gluon plasma" ) [l) . A multitude of relevant physics observables provides 
indication {2-8] that this phase transition might first occur in the energy range of the CERN 
SPS (6.3 GeV < y/s NN < 17.3 GeV): the domain of the present study. 

While the above evidence is based on inclusive observables, additional insight can be 
obtained by looking at event-by-event fluctuations. Enhanced fluctuations are a general 
feature of phase transitions. Recent lattice QCD calculations report indications of a QCD 
critical point at finite baryo-chemical potential, reflected in a steep rise of quark number 
density fluctuation (quark "susceptibility") |2H|. In addition, the changed correlation pat- 
terns in particle production, expected where new degrees of freedom emerge, are accessible 
via fluctuation observables 1121 . 



Numerous such observables have been explored previously [12l-ll4|. in part with inconclu- 
sive outcome. This resulted from uncontrolled sources of background fluctuations, notably 
variation of the reaction volume due to concurrent impact geometry fluctuations, finite num- 
ber statistics, or from an obliteration of the signals during the hadron-resonance expansion 
phase that follows hadron formation at the parton-hadron phase boundary (notably so for 



electric charge correlations 
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The present study of hadron ratio fluctuations overcomes the mentioned difficulties to 
a certain extent. Effects from volume fluctuations are minimized, as a particle ratio is 
an intensive quantity and the average hadron ratios change only marginally in the analysed 
centrality interval. Fluctuations from finite number statistics, as well as effects of the limited 
detector resolution in the particle identification (PID) method are removed by subtracting 
a mixed event reference. Thus, dynamical fluctuations are studied. 

For kaon-to-proton number ratio fluctuations, the conserved charges strangeness and 
baryon number are carried by heavy particles, making their dispersion in momentum space 
smaller compared to the case of electric charge fluctuations, where the Q-value of resonance 
decays creates a noticeable difference between the original charge distribution and that of 
the finally observed pions [13|. No such transport process (as, e.g. a resonance feeding into 
K + + p) is known in the case discussed here. 

The results presented here are the first attempt to experimentally probe the baryon 
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number-strangeness correlation, which undergoes a massive change at the deconfinement 
phase transition. Above that transition, where strange quarks (S = — 1, -6 = 1/3) are the 
relevant degrees of freedom, strangeness S can only exist in direct conjunction with baryon 
number B. In a hadron gas, kaons (S = —1, B = 0) allow for strangeness number unrelated 
to baryon number. An appropriate correlation coefficient Cbs 15[ has been proposed which 
is related to the above mentioned quark number density susceptibililties at finite baryo- 
chemical potential ^b recently predicted by lattice QCD [9-ll|. The latter exhibit steep 
maxima at T = 150 MeV and //b ~ 400 MeV, conditions that are reached in A+A collisions 
in the energy range under investigation in this analysis. An important contribution to this 
coefficient are the kaon-proton correlations, that are experimentally probed for the first time 
in the present data analysis. In a dense hadronic medium, in contrast to a gas of uncorrelated 
hadrons, a variety of processes may lead to a correlation of baryon number and strangeness, 
e.g. the associated production of hyperons and strange mesons. The effect of these hadronic 
correlations on Cbs was studied in the hadronic transport model UrQMD 2lH23| and found 
to be small and constant in the SPS energy range. 

The following analysis is based on data from Pb+Pb collisions recorded at the CERN 
SPS at five energies, y/s NN = 6.3, 7.6, 8.7, 12.3 and 17.3 GeV. The 3.5% most central events 
were selected, also avoiding an influence of the change in the inclusive (K)/(p) ratio with 
centrality. 

The NA49 detector [16] uses 4 large volume time projection chambers (TPC) for tracking 
and particle identification (PID) via their specific energy loss (dE/dx) in the TPC gas. 
During the SPS energy scan program, care was taken to keep the acceptance approximately 
constant with respect to midrapidity by setting the magnetic field strength proportional to 
the beam momentum. To ensure best resolution and stability of the dE/dx measurement, 
only tracks within the acceptance of the large main TPCs were accepted for the analysis. 
Further quality cuts are applied on the distance of closest approach of the extrapolated 
particle trajectory back to the main vertex and on the number of measured dE/dx samples 
on the track. The resulting acceptance as a function of center-of-mass rapidity normalized by 
beam rapidity, and transverse momentum p^< is depicted in Fig. [TJ The azimuthal acceptance 
is described in [l7J and the comprehensive acceptance tables to be used in simulations can be 



found in [18(. At each energy, the largest possible acceptance was used to ensure the highest 
statistical significance. Average uncorrected multiplicities within this acceptance are given 
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FIG. 1: The acceptance for kaons and protons, as used in the present analysis, as a function 
of transverse momentum and the center-of-mass rapidity, normalized by the corresponding beam 
rapidity, for all analysed energies. Lines delimit the regions in which particles can be identified. 
Limits are due to geometric acceptance and available statistics, the latter predominantly at large 
momenta p and transverse momenta p^. 

m table E 

The analysis procedure is similar to the one employed by NA49 in the study of fluctuations 



of the K/ir and p/n ratios 171 ] . The dE/dx resolution of approximately 4% allows for a 
statistical determination of the inclusive particle yields through a \ 2 fit to the dE/dx spectra 
in momentum space bins. Probability density functions f m for dE/dx and momentum 
F m from this inclusive analysis are then used as input for the event-wise particle ratio 
determination. The dE/dx distributions overlap for different hadron species m, making a 
simple particle counting impossible. Therefore, an unbinned likelihood method as introduced 
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v^W(GeV) <K+) 00 (P) <P> 

6.3 5 1 28 

7.6 8 2 35 

8.7 12 4 41 
12.3 22 9 53 1 
17.3 34 20 72 3 



TABLE I: Average uncorrected multiplicities within the acceptance of the present study. 



in 



19 



20J is used. 

In each event, the likelihood function L is obtained by multiplying the probabilities of the 
n particles in the event ((n) ~ 60 at y/s NN = 6.3 GeV and (n) « 600 at y/s NN = 17.3 GeV): 



n 

i=i 



(pi) fm {fh {dE/dx)i 



(1) 



where 9 m are the relative abundances, constructed such that J2 m ©m = 1- They are the only 
free parameters in the subsequent likelihood maximization. The event-wise hadron ratios 
are then calculated from the fitted values of m . 

The fluctuations of the event-wise K/p ratio are quantified using the scaled dispersion 



a := J Var(K/p)/(K/p). The measured value in the data is denoted Odata- Figures |2] and [3] 
show the event-by-event distribution of the (K + + K~)/(p + p) ratio for y/s NN = 17.3 and 
6.3 GeV. To quantify finite-number statistics and PID resolution effects, a reference sample of 
mixed events is constructed. This sample is made such as to preserve the original multiplicity 
distribution, with no two tracks in a mixed event taken from the same physical event. As 
described in detail in [l^] , the mixed events contain no correlation due to physical processes, 
but effects from finite number statistics remain. The measured dE/dx information is still 
attached to the individual particles, so that the likelihood method can be applied to the 
mixed events in the same way as to the original events. Thus, the effect of the dE/dx 
resolution on the extracted particle ratios is reproduced by the mixed events. The event- 
by-event distribution from mixed events is also displayed in Figs. [2] and [31 Their scaled 
dispersion is denoted cr mix . 



7 



1.2 



x 
E 

T3 



(2 0.8 



H — I — | — I — I — I — I — | — I — I — I — I — | — 1 n l MM 




0.5 1 1.5 

Event-by-event (K + +K")/(p+p) 



FIG. 2: (Color online) Event-by-event distribution of the (K + + K~)/(p + p) ratio at ^s NN = 
17.3 GeV for real data events compared to the mixed event reference. The lower panel shows the 
ratio data/mixed, where the convex shape indicates negative dynamical fluctuations. Only in the 
ratio plot, for better readability, statistically insignificant bins are not shown. 



The dynamical fluctuations can now be constructed as the quadratic difference 17l.ll9l.l20l]: 

sign (■ 



0"d 



yn 



data 



data 



(2) 



At y/s NN = 17.3 GeV (Fig. [2]) the data show a narrower distribution compared to mixed 
events, corresponding to a negative Od yn - This is further visualized in the ratio between data 
and mixed event distributions shown in the lower panel of Fig. [2] . Conversely, at 6.3 GeV 
(Fig. [3]) o"data > °mix an d thus o"d y n > 0. Following from the less symmetrical shape of the 
event-by-event distribution here, this result is not so straight-forwardly inferred from the 
ratio plot. 

The method described above has been successfully used and thoroughly tested in the anal- 
ysis of K/tt and p/n fluctuations [17|, and the same extensive quality checks were applied in 
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FIG. 3: (Color online) Event-by-event distribution of the (K + + K~)/(p + p) ratio at ^/s NN = 
6.3 GeV for real data events compared to the mixed event reference. The lower panel shows the 
ratio data/mixed, where the concave shape indicates positive dynamical fluctuations. 

the present analysis. For instance, outlying events with very small or high K/p ratios were 
found to contribute only modestly to the reported fluctuation signal: Consistent with jl?! ]. 
the signal changes by less than 1% when removing the high or low tails of the K/p distri- 
butions that correspond to 1% of the events. Special care was taken in the present analysis 
to check whether a correlation remains from the fit when extracting the K/p ratio, as the 



dE/dx based kaon-proton separation is smaller than for the kaon-pion case. 



an influence, events generated in the hadronic transport model UrQMD 21 



b exclude such 



231 ] were used 



to study the effect of the dE/dx resolution. The results of a direct model calculation are 
in agreement with those obtained using in addition a simulation of the NA49 TPC dE/dx 
response and employing the fit procedure described above. The difference between the two 
methods amounts to 1.5% at most and is taken into account in the systematic error. 
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FIG. 4: (Color online) Energy dependence of the event-by-event dynamical fluctuations of the 
(K + +K _ )/(p+p) and the K + /p ratios. Symbols represent the NA49 measurements with statistical 
and systematic (braces) uncertainties. Calculations within the UrQMD and HSD transport models, 
processed through an NA49 acceptance filter are represented by lines. The statistical errors on the 
model results decrease from approximately 1.5% at 6.3 GeV to 0.5% at 17.3 GeV. 



The stability of the results was tested under variation o 



the track selection criteria, 



acceptance and event selection, as described in detail in [l7|- The changes due to these 
variations are also represented in the systematic errors. The results proved to be stable 
under small variations of the chosen acceptance. However, the acceptance tables provided 
in 18j should be used for model comparisons. 

(Tdyn has been evaluated for two ratios: (K + + K~) /(p + p) and K + /p. The former probes 
correlations between all four involved hadron species. Only K + /p can be studied among the 
separate charge ratios, as the K~ and p yields (see table HI) vanish at low energies. For the 
same reason, the two observables are expected to converge at low energies. The excitation 
functions of <J dyn &re shown in Fig. HI For both cases, the dynamical ratio fluctuations 
change from positive (enhanced fluctuations compared to mixed events) at low SPS energies 
to negative values at the higher SPS energies. At high and low energies, <Td yn agrees for the 
two studied ratios, but disagrees at y/s NN = 7.6 and 8.7 GeV. Recent preliminary results 
from the STAR collaboration j^J also find negative values and indicate a weak energy 
dependence between y/s NN = 17.3 and 200 GeV. 
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23| and HSD [23, Um are 



Calculations in the hadronic transport models UrQMD 21 
shown for comparison. The NA49 experimental acceptance was used in the model studies. 
(Xdyn shows only a weak energy dependence in both models. The striking change of sign and 
the strong energy dependence seen in the data are not reproduced. In UrQMD, both charge 
combinations were evaluated and have a constant difference over the studied energy range, 
another feature in contrast to the data. 

Unlike the previous results on K/tt and (p + p) /n fluctuations, for which attempts of 
an explanation exist, the present data can not easily be understood. We recall that the 
(p + p) /n fluctuations were explained in hadronic models as a result of the proton-pion 
correlation due to resonance decay 



jl7 . 127I and that the rise in cTd yn for K/ir was suggested 



to be due to scaling properties of the observable <7d yn itself 



to the onset of deconfinement 



or might even be connected 



29]. 



As no known resonance feeds into positively charged kaons and protons, another source 
of correlation has to change at the energy where a^ yTl switches sign, deviating from the 
transport model calculations. The change of sign is also incompatible with all scalings 
suggested in 2i| and thus indicative of a change in the underlying correlation physics. 

In the baryon-strangeness correlation, a rapid change is expected at the deconfinement 



phase transition 30_|, |31( for which indications were found in the same energy region in rapid 
changes of several hadron production properties 5KZ|. The observed energy dependence 
qualitatively supports the scenario of a change in the baryon-strangeness correlation, but 
the exact contribution to Cbs of crd yn ((K + + K~)/(p + p)) and <7d yn (K + /p), respectively is 
still under discussion [32 ] . 

In summary, we present a first measurement of the dynamical fluctuations of the kaon- 
to-proton number ratio at the SPS energies. Both (K + +K _ )/(p + p) and K + /p fluctuations 
show a non-trivial excitation function that is not reproduced in the hadronic models UrQMD 
and HSD, potentially pointing to a change in the baryon number-strangeness correlation at 
y/s NN ~ 8 GeV. Although a connection between the kaon-to-proton ratio and Cbs seems 
reasonable, and the latter is suggested as a unique test for the basic degrees of freedom 
in the probed matter, the detailed connection between our measurement and Cbs an d its 
interpretation require further theoretical studies. 

This work was supported by the US Department of Energy Grant DE-FG03- 
97ER41020/A000, the Bundesministerium fur Bildung und Forschung, Germany, German 



11 



Research Foundation (grants GA 1480/2-1 and STO 205/4-1), the Polish Ministry of Sci- 
ence and Higher Education (1 P03B 006 30, 1 P03B 127 30, 0297/B/H03/2007/33, N N202 
078735, N N202 204638), the Hungarian Scientific Research Foundation (T068506), the Bul- 
garian National Science Fund (Ph-09/05), the Croatian Ministry of Science, Education and 
Sport (Project 098-0982887-2878), and Stichting FOM, the Netherlands. 



F. Karsch and E. Laermann, arXiv:hep-lat/0305025 



R. Arnaldi et al. [NA60 Collaboration], Eur. Phys. J. C 61, 711 (2009) [arXiv:0812.3053l 
[nucl-ex]]. 



U. W. Heinz and M. Jacob, arXiv:nucl-th/0002042 



M. Gazdzicki and M. I. Gorenstein, Acta Phys. Polon. B 30, 2705 (1999) 
[arXiv:hep-ph/9803462] . 

S. V. Afanasiev et al. [The NA49 Collaboration], Phys. Rev. C 66, 054902 (2002) 
[arXiv:nucl-ex/0205002] . 

C. Alt et al. [NA49 Collaboration], Phys. Rev. C 77, 024903 (2008) [arXiv:0710.0TT8l [nucl-ex]] . 
V. Friese, PoS CPOD09, 005 (2009) 



M. Gazdzicki et al. [NA49 Collaboration], J. Phys. G 30, S701 (2004) |arXiv:nucl-ex/0403023 1. 

Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004) [arXiv:hep-lat /0402006] . 

S. Ejiri, C. R. Allton, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann and C. Schmidt, 

Prog. Theor. Phys. Suppl. 153, 118 (2004) |arXiv:hep-lat/0312006| . 

F. Karsch, PoS CPOD07, 026 (2007) |arXiv:0711.0656l [hepdat]]. 

V. Koch, arXiv:0810.2520 [nucl-th]. 

J. Zaranek, Phys. Rev. C 66, 024905 (2002) |arXiv:hep-ph/0111228| . 

C. Alt et al. [NA49 Collaboration], Phys. Rev. C 70, 064903 (2004) l arXiv:nucl-ex/0406013| . 
V. Koch, A. Majumder and J. Randrup, Phys. Rev. Lett. 95, 182301 (2005) 
[arXiv:nucl-th/0505052| . 

S. Afanasiev et al. [NA49 Collaboration], Nucl. Instrum. Meth. A 430, 210 (1999). 

C. Alt et al. [NA49 Collaboration], Phys. Rev. C 79, 044910 (2009) [arXiv: 0808 .12371 [nucl-ex]] . 

https: / / edms.cern.ch /document / 984431 / 1 

S. V. Afanasiev et al. [NA49 Collaboration], Phys. Rev. Lett. 86, 1965 (2001) 



12 



[arXiv:hep-ex /0009053| . 
[20] M. Gazdzicki, Nucl. Instrum. Meth. A 345, 148 (1994). 
[21] M. Bleicher et al, J. Phys. G 25, 1859 (1999) |arXiv:hep-ph/9909407] . 
[22] S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998) [arXiv:nucl-th/9803035] . 
[23] H. Petersen, M. Bleicher, S. A. Bass and H. Stocker, arXiv:0805.0567l [hep-ph]. 
[24] J. Tian [STAR Collaboration], J. Phys. G 37, 094044 (2010). 
[25] W. Ehehalt, W. Cassing, Nucl. Phys. A602, 449-486 (1996). 

[26] V. P. Konchakovski, M. Hauer, M. I. Gorenstein, E. L. Bratkovskaya, J. Phys. G 36, 125106 

(2009). |arXiv:0906.3229l [nucl-th]]. 
[27] T. Anticic et al [NA49 Collaboration], PoS CPOD2009, 029 (2009) |arXiv:0910.0558l [nucl- 

ex]]. 

[28] V. Koch and T. Schuster, Phys. Rev. C 81, 034910 (2010) |arXiv:0911.ll60l [nucl-th]]. 

[29] M. I. Gorenstein, M. Gazdzicki and O. S. Zozulya, Phys. Lett. B 585, 237 (2004) 

|arXiv:hep-ph/0309142] . 
[30] M. Cheng et al, Phys. Rev. D 79, 074505 (2009) |arXiv:0811.1006l [hep-lat]]. 
[31] C. Schmidt, PoS CPOD2009, 024 (2009) |arXiv:0910.432T1 [hep-lat]]. 
[32] V. Koch, private communication. 



13 



